Holographic Dark Energy with Curvature 



O 

(N 

CU 

<D 

(N 

o 

u 

43 

Oh- 
i 

o 

CO 

03 



> 
O 

On 

in 

o\ 
o 
o 



13 



Victor H. Cardena:£] and Roberto G. Perej]] 
Departamento de Fisica y Astronomia, Universidad de Valparaiso, Gran Bretana 1111, Valparaiso, Chile 

In this paper we consider an holographic model of dark energy, where the length scale is the 
Hubble radius, in a non flat geometry. The model contains the possibility to alleviate the cosmic 
coincidence problem, and also incorporate a mechanism to obtain the transition from decelerated to 
an accelerated expansion regime. We derive an analytic form for the Hubble parameter in a non flat 
universe, and using it, we perform a Bayesian analysis of this model using SNIa, BAO and CMB 
data. We find from this analysis that the data favored a small value for Qk, however high enough 
to still produce cosmological consequences. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 

One of the most important problems in cosmology is 
to explain the fact that the universe today is in a phase 
of accelerated expansion. In 1998, two teams studying 
distant Type la supernovae presented independent evi- 
dence that the expansion of the Universe is speeding up 
[H,i|. The physical origin of cosmic acceleration remains 
a deep mystery. Since the expansion is speeding up, we 
are faced with two possibilities, either of which would 
have profound implications for our understanding of the 
cosmos and the laws of physics. The first is that 75 % 
of the energy density of the Universe exists in a new 
form with large negative pressure, called dark energy. 
The other possibility is that General Relativity breaks 
down on cosmological scales and must be replaced with 
a more complete theory of gravity. Among the models 
proposed so far for dark energy, the Holographic Dark 
Energy (HDE) scenario emerge as one of the few that 
has a deep physical substratum. 

The idea of Holography was first suggested by Beken- 
stein Q, by generalizing the results obtained from black 
hole physics, and postulating that the maximum entropy 
in a box of volume L 3 should growth as the area of the 
box. In this context, 'tHooft Q and Susskind @ stressed 
that this principle has to be considered as a recipe for a 
quantum gravity theory. The idea soon finds its way 
in cosmology The Holographic Principle can be de- 
scribed as follows: if we want to reconcile quantum me- 
chanics with gravity, we have to assume that the observ- 
able degrees of freedom of the universe are projections 
coming from a two-dimensional surface, where the infor- 
mation is stored. 

Considering that particle physics can be successfully 
described by an effective field theory with an ultraviolet 
(UV) cutoff, and black hole physics impose an infrared 
(IR) cutoff to the energy density involved, Cohen et al. Q 
suggested that the dark energy should satisfy this prin- 



ciple, proposing a relationship between cutoff as 

L 3 P a < LMp , (1) 

where p\ is the dark energy density, L the size of the 
system, and M p — (8-kG)^ 1 ^ 2 is the reduced Planck mass. 
In [8J the author notice that using L as the size of the 
observed universe (Ltoday — 10 Gy) yields a value p\ ~ 
1CU 10 eV 4 , which would explain the observed value of 
the dark energy density. However, Hsu [9f demonstrated 
that although identifying L with the Hubble radius H~ x 
leads to a correct order of magnitude for the dark energy 
density, written p\ — 3c 2 H 2 M 2 saturating (fTJ) with c 
constant, this leads to a wrong equation of state. Based 
on this, Li [loj proposed instead to identify L with the 
future event horizon. 

Soon after, Pavon and Zimdahl llj showed that is still 
possible to consider the Hubble radius as a len gth scale in 
the context of interacting dark energy models [12| , which 
allows, as a bonus, to alleviate the cosmic coincidence 
problem: why the energy density of matter and dark en- 
ergy are of the same order of magnitude just today. This 
time however, pa = 3c 2 H 2 M 2 can actually describe the 
dark energy density assuming that c is a time dependent 
function. The analysis performed in that paper was ex- 
tended in [l3[ for non-flat cosmologies, focusing in the 
evolution of the ratio of the energy densities r = p m / PA- 
An statistical analysis of this model in a flat universe was 
performed recently in 14 1. 
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In this paper we consider an interacting holographic 
dark energy model, using the Hubble radius as the length 
scale, in a universe with arbitrary curvature. The model 
not only contains the possibility to alleviate the coin- 
cidence problem, but also incorporate a mechanism to 
obtain the transition from decelerated to an accelerated 
expansion regime. We are primary interested in study the 
effects of the inclusion of a small curvature in the theo- 
retical parameters in light of the observational data. We 
derive an analytic form for the Hubble parameter in a non 
flat universe, and using it, we perform a Bayesian analy- 
sis of this model using supernovae data, Baryon Acoustic 
Oscillation (BAO) distance, and the CMB shift param- 
eter. We find from this analysis that the data favored 
a small value for f2fe, however high enough to produce 
important cosmological consequences [lj| . 
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In the next section we review the model clarifying the 
assumptions made in [13| and [14]]. In section III we 
derive the flat universe results, some of them well known, 
and others new. The main contribution is explained in 
section IV where is derived the general solution for the 
non flat case. Next, we perform the bayesian analysis 
using observational data, and we end the paper with some 
comments. 



II. AN EFFECTIVE EQUATION OF STATE FOR 
DARK ENERGY 

The Einstein field equations for a homogeneous and 
isotropic universe are the Friedman equation 



We can also rewrite this equation by using ([5]) and ([6]) 
obtaining 



// _ I d In p A 

independent of the value of w\. 



(10) 



III. THE HOLOGRAPHIC MODEL FOR A FLAT 
UNIVERSE 

As a warm up and also to introduce the formalism, let 
us discuss the flat model first. This section is a critical 
review of previous works 0] > 0] y 0] ■ 



H 2 + ^ = 



1 



3M 2 



P: 



(2) 



where p is the energy density, H = a/a and a(t) is the 
scale factor, k = +1,0,-1 for a closed, fiat, and open 
geometries respectively, and 



H 



1 



2M 2 



(p + P) + -, 



(3) 



where P is the pressure, and where a dot means cosmic 
time derivative. The model emerges by assuming that 
the energy density p, which satisfy the usual energy con- 
servation equation 



p + 3H(p + P) =0, 



(4) 



can be separated into two components; dark matter and 
dark energy, so p = p m + p\ and P = p m +p\. Assuming 
explicitly the equation of state for dark matter, p rn = 0, 
and for dark energy p\ = w\Pa, and further, considering 
that neither component conserves separately, we obtain 



Pa + 3Hp A (l + w\) 
Pm + 3Hp rn 



-Q, 



(5) 
(6) 



where Q is an arbitrary interaction term between dark 
matter and dark energy. From ([5]) we can take the Q 
term to the left hand side and write 



which looks like the energy conservation equation for 
dark energy but with an effective pressure given by 



eff 

Pa =Pa 



Q_ 

3H' 



(8) 



Then, we can write an effective equation of state for this 
dark energy as 



ui 



eff 

•Jf _ Pa 



WA 



'SHpA 



(9) 



A. Constant c 

As we stated in the introduction, the corresponding 
vacuum energy density can be written from the holo- 
graphic dark energy as 



Pa = 3c 2 M 2 H 2 , 



(11) 



where c is considered constant. Combining this with the 
Friedman equation enable us to write the energy density 
of dark matter 



Pm = 3M 2 H 2 {1 



(12) 



To protect the positivity of the dark matter energy den- 
sity, c < 1 is required. Clearly, if c is constant and the 
two components evolve independently (i.e. interaction 
Q = 0), then p m ~ H 2 ~ a -3 as the usual dark mat- 
ter, but this implies also that pa — a~ 3 , an evolution 
that can not generate accelerated expansion. This is the 
result found by Hsu [|, 0. 

However, as first stated in 01, if we consider an inter- 
action between dark matter and dark energy (see Eqs.© 
and ([6])), we can obtain an equation of state for dark en- 
ergy different from that of dust. This is exactly what we 
have derived in the previous section in Eq. (|10[) . Because 
H = a/a and using Eq.Q, we can find the dark energy 
density as a function of the scale factor 



PA 



1-c 2 



-3(l-c 2 ) 



(7) In this case the deceleration parameter is 



q = \{l-ic 2 ). 



(13) 



(14) 



To obtain an accelerated expansion universe, i.e. q < 
0, and to protect the positivity of dark matter energy 
density, one obtains a range of values the constant c can 
take: 1/3 < c 2 < 1. The effective equation of state of 
vacuum energy density is 



eff _ 



_ 1 d In pa 
3 din a 

-c 2 . 



(15) 
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Relation (|T3|) indicates that the ratio r = p m /p\ = 
(1 — c 2 )/c 2 is a constant. Although it is desirable to 
get a nearly constant r parameter, solving the coinci- 
dence problem, getting a constant value, produces an- 
other problem: there is no transition from decelerated to 
an accelerated expansion regime. 



B. Variable c 

As we notice at the end of the previous section it is 
clear that when c is constant, no transition from deceler- 
ated expansion to accelerated expansion can be realized, 
a transition observationally tested. A possible way to 
overcome this difficulty was proposed in (llj . The idea 
is to consider c as a time dependent parameter, enabling 
also that r = (1 — c 2 )/c 2 to be. Since r is expected to de- 
crease with time, this implies that c 2 increases, keeping 
the relation c 2 (l + r) = 1. 

In this case 



p A = 3c 2 {t)M 2 P H 2 . 



(16) 



and also, using the Friedman equation, we can write the 
relation 



p m = ZM 2 P {l-c 2 (t))H 2 . 



(17) 



Notice that, although the redshift dependence of p\ is 
not proportional to H 2 , no further constraint on c 2 (z) 
is derived from the Friedman equation, because p m in 
(fl~7|) is defined to satisfy it automatically. Again, it is to 
protect the positivity of the energy density of cold dark 
matter that c(t) 2 < 1 is required. From the conservation 
equation of cold dark matter Eq. and the Friedmann 
equation, one can write in a flat universe 



(1 



^ln-ffO) 

dz 



w A c 2 {z)). 



(18) 



To solve Eq. (fT5|) . one has to assume a concrete form of 
the function c(z). After a simple calculation, one can 
also write the deceleration parameter 



;(l-3c 2 (z)). 



(19) 



a relationship first written in [11(, that shows how the 
equation of state parameter depends explicitly on the 
interaction and the c 2 (t) function. 

If we write Q — Tp\, with T being some particle pro- 
duction rate, we can write a solution for ([6]) 



Pn 



Pi 

V a ■ 







) exp 









(22) 



where the subscript indicates today. An explicit model 
for a flat universe and interaction was proposed in f]~3j j . 
where it was suggested that 



r 

7h 



3/3 



(23) 



in such a way that (f2"2"j) can be solved exactly. In this 
case 



Pv 

p° 



(a y 

\T) exp 



1 



(24) 



If we assumed that c 2 is constant, then pa evolve in the 
same way as p m , so we can explicitly write down the 
Hubble parameter from the Friedman equation. The de- 
celeration parameter can also be computed being 



y 2 T \ao 



(25) 



An identical expression for q was derived in 14|, where 
instead a explicit form of c 2 {z) was given 



(l + z) a 



(26) 



The equivalence between the models can be easily ex- 
plained using the equation (|21[) . The analysis performed 
in [14| consider a flat universe k = with interac- 
tion and is explicitly used that w\ = — 1. This last 
assumption implies a relation between the interaction 
Q and c 2 . Actually, up to a constant, the relation is 
c 2 = Q/(3rH pa). The analysis performed both in [l3[ 
and [14j correspond to this model, considering interaction 
between the dark constituents in a flat universe, keeping 
the value c 2 (l + r) = 1. 



One finds that once < c 2 (z) < 1 is time dependent, 
possible transition from deceleration expansion to accel- 
erated expansion can be realized. The analogous expres- 
sion for the effective equation of state in this case, can 
be derived using Eqs. ([2]) and (|16p . obtaining 



w 



eff 



W A C 



ldlnc 2 



3 din a 

Combining this equation with ^ we obtain 



1 + r 



3H P a 3Hc 2 



(20) 



(21) 



IV. NON FLAT MODEL 

In this section we proceed to generalize the previous 
work to non flat universes. Although it is often said that 
the inflationary scenario predicts a flat universe [l6| , the 
situation has not been settled down yet, partly because 
the observations are not conclusive and because in the- 
ory rife = (p,k = —k/a 2 H 2 in this paper) is not the 
only solution [17| . Moreover, in the context of parameter 
estimations, considering the curvature as a new free pa- 
rameter, leads to an increment in the uncertainties in all 
the others cosmological parameters. In particular, in the 
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context of models with dynamical dark energy, a number 
of works have discussed the degeneracies between w(z) 
and flk [HI, [HI; (Ifj- Therefore, the precision in obser- 
vational cosmology requires the inclusion of the curva- 
ture, even if its contribution is small, in any theoretical 
study of dark energy. In this section c 2 is a function of 
time and fly ^ 0. 

Let us start with Eq.Q. Replacing the definitions of 
both p m and pa we found 



H 



(27) 



Using that a = (1 + z) 1 we can rewrite this equation 
using the redshift as the independent variable 

(1 + z) dH - OTj2n 1 -■ - 2 



dz 

which is the generalization of (|T8 
pTjl follows in the same way as 



3iT z (l + WAC ^ + k(l + z ) 2 1 (28) 
The generalization of 



U>A 



1 



fli 



Q 



3H PA 



3 



3Hc 2 



(29) 



that contains the obvious limit Slfc = result (|2"Tj) . and 
also coincides with the result given in (l3| (see equation 
(26) in that paper). To show the equivalence we use the 
identities fl m + flk + flA — 1 and r — fl m /fl\. The analog 
of equation (j2TJ)) in this case shows how the presence of 
curvature changes the dynamics. The effective equation 
of state parameter when w\ = — 1 takes the form 



w 



eff 



flk 1 d In c 2 
3 3 din a 



(30) 



The presence of curvature can ameliorate the negativity 
of w e ^ or, if it is positive, it can be even more negative, 
crossing the phantom line. However, the statistical anal- 
ysis performed on this model, that we describe below, 
shows that because the curvature parameter turns out to 
be very small (see results in Table I), the actual value 
of flk has a little impact on the dynamics of iv^f (see 
Figured]). 

An explicit solution of (f2"5f is possible when we fix the 
bare equation of state parameter wa = — 1 



TABLE I: The best fit values for the free parameters in the 
case of a non flat holographic dark universe, a and (3 charac- 
terize the function c (z) defined in (|26l) . and flk is the curva- 



ture parameter. 


We specify the data 


set used in 


each . 


Data Set 


Xmin a ± 0.02 


/3L0.08 


fl k ± 0.009 


SN 


310.866 1.089 


1.513 


-0.5873 


SN+BAO 


310.406 0.732 


3.161 


0.3876 


SN+BAO+CMB 312.848 0.728 


1.219 


-0.0449 



Having the Hubble parameter in terms of the redshift, 
we can compute the luminosity distance and confront the 
model with observations of supernovae (SN) [22J. In this 
work, we use the Supernova Cosmology Project Union 
sample [2l|, having 307 SN distributed over the range 
0.015 < z < 1.551. We fit the (theoretical) distance 
modulus fi(z)th defined by 



K z )th = 



H d L (z) 



Mo, 



(33) 



to the observational ones ji{z) b s - Here Hq = 100/ikm 
s _1 Mpc -1 is the Hubble constant and the luminosity 



distance is defined by d^(z) — (1 



r{z) 



H i/\fJ, 



-Smn^f\fTk~ 



z)r(z) where 
dz' 



o H(z'Y 



(34) 



and /io = 42.38 — 5\og 10 h. Sinn(a:) = sin a;, x, sinh x for 
flk < 0, flk =0, and fl k > respectively. Considering 
a, (3 and flk as free parameters, we get a Xmin = 310.866 
with the best fit values displayed in Table I. 

The second major input for parameter determination 
comes from the baryon acoustic oscillations (BAO) de- 
tected by Eisenstein et al. [23| ■ In our work, we add the 
following term to the % 2 of the model: 



Xbao 



(A - A BAO ) 



where A is a distance parameter defined by 

-l 1/3 



A = 



\/fl m H 2 

czbao 



r (zbao) 



czbao 
H(zbao) 



(35) 



(36) 



3/3 



30 flk 

6 a+ Tv> 



z?e 



H\z) = H 2 (l 



a a a a 



x (31) 



where L(a,z) is the incomplete gamma function or plica 
function defined as 



and A BAO = 0.469, a A = 0.017, and z B ao = 0.35. 
Considering both supernovae and BAO data, we obtain 
Xmin — 310.406 with the best fit values displayed in the 
second row of Table I. 

The CMB shift parameter R is given by [H[ 



R(z. t 



(37) 



T(a,z) 



t a - x e- l dt. 



(32) 



Expression (l3~Tj) for fly. = reduces to the expression first 
written in [l3| and contrasted with observations in [i"4| . 



Here the redshift z* (the decoupling epoch of photons) is 
obtained by using the fitting function [25j 



2* = 1048 [l + 0.00124(fi h /i 2 



-°- 73s ] [i+ 91 (fi m h 2 r], 

(38) 
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FIG. 1: These are the confidence contours in the plane a — j3 
with la and 2a regions. The central dot shows the best fit 
parameter for the joint analysis using SN + BAO + CMB 
indicated in the third line of table I. This case correspond to 
a curvature parameter of fifc = —0.0449 ± 0.009. 



where the functions g\ and 32 are given as 



5 i = 0.0783(f^ 2 )- - 238 (1 + 39.5(ft b /i 2 ) a763 ) (39) 



92 



0.560 (1 + 21. l(Q b h 



(40) 



The WMAP-7 year CMB data alone yields R(z*) = 
1.726 ± 0.018 26]. Defining the corresponding Xcmb 
as 



Xcmb 



(fl(z,) - 1.726) ; 
0.018 2 



(41) 



one can deduce constraints on a, ft and f2fc. A joint anal- 
ysis using SN+BAO+CMB leads to the best fit values 
showed in Table I. As is well known, we notice the effect of 
incorporating the BAO and CMB data to the supernova, 
leading to a decreasing value for Qj.. As in the ACDM 
model, when we use only the SNIa data, in particular the 
Union compilation, the fit sug gest a huge curvature con- 
tribution (see figure 11 in [2l[ and the discussion). Only 
when we incorporate BAO and CMB constraint, the cur- 
vature parameter tends to smaller values. As is also well 



known, the strongest observational indication for a fiat 
universe is the CMB data. A gaussian approximation to 
the confidence contour is also shown in Figure [TJ where 
we show the plane a — j3 for = —0.0449. Using the 
best fit parameters we show in Figure [2] a plot of the 
resulting effective equation of state ([30]) . 

By choosing a cosmological constant equation of state 
w\ = —1, we get a model that, in agreement with ob- 
servations, seems to cross the phantom line, leading to 



T for z < 1. 



-1.2 



FIG. 2: Using the best fit parameters found in the analysis 
(see Table I), we plot the effective equation of state, expression 
(|20[) . for our model. Notice that even if the bare equation of 
state parameter iua — — 1, the combined data is consistent 
with an effective equation of state that cross the cosmological 
constant line towards phantom like models. 



V. DISCUSSION 

We consider a non flat universe version of the holo- 
graphic dark energy model studied first in 13 1 and [Tij ]. 
in which the interaction between dark matter and dark 
energy enable us to describe the transition from a de- 
celerated to an accelerated expansion regime, together 
with alleviating the cosmic coincidence problem. Taking 
into account the spatial curvature, we were able to find 
an analytic expression for the Hubble parameter, which 
was used to perform a direct confrontation to observa- 
tional data as SNIa, BAO and CMB shift parameter. 
This analysis shows that our model is compatible with 
current precision data. The best fit value for the curva- 
ture obtained, although small, it is well within the order 
where its presence can modify for example, the number 
of e-folding of inflation which changes dramatically the 
early stages of inflation [15[ . 
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